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lonospheric Disturbances Produced by 


Severe Thunderstorms 


KENNETH DAVIES AND JOHN E. JONES 
Space Environment Laboratory, Environmental Research Laboratories, 


Boulder, Colo. 


ABSTRACT. Radio observations near Oklahoma City during the spring of 1970 
have shown ionospheric effects that occur within a radius of 200 km from severe 
thunderstorms having cloud tops above 12 km (40,000 ft). These effects frequently 
take the form of acoustic-gravity waves that appear to originate in the thunder- 
storm cells and propagate to F-region heights. The power in the waves usually 
has peaks at periods near 3.5 and 4.5 min, although periods between 7 and 10 min 
have been observed occasionally. Three-dimensional velocity vectors have been 
calculated for several ionospheric disturbances; and by means of acoustic wave 
ray-tracing, particular groups of thunderstorm cells have been identified as their 
sources. For the 7- to 10-min waves, an apparent inconsistency between the period 
and the velocity gives evidence of neutral winds with speeds of at least 70 m/s in 
the F region. We estimate that the power required to produce waves of the type 
observed is about 50 MW, one-millionth of the power generated by a large storm. 


1. INTRODUCTION 


We discuss certain types of wavelike disturbances 
observed on ionospheric radio echoes. We show that 
these disturbances are closely associated with the 
occurrence of severe thunderstorms within a hori- 
zontal distance of about 200 km from the ionospheric 
reflection point. The evidence suggests that the 
observed phenomena result from the interaction 
between ionospheric charged particles and the 
motion of neutral molecules caused by the passage 
of acoustic type waves with periods between 3 and 
5 min. 

In section 2, we discuss the evidence in the 
literature for ionospheric effects thought to be asso- 
ciated with tropospheric weather. Section 3 contains 
some basic radio propagation theory required for 
understanding the probing technique with spaced 
transmitters and also contains the description of 
the method of data analysis. Some physical theory of 
the propagation of acoustic gravity waves is pre- 


sented in section 4 together with the results of ray 
tracing in the 1962 U.S. Standard Atmosphere 
[COESA(U.S. Committee on Extension to the Stand- 
ard Atmosphere) 1962] with and without horizontal 
winds. In section 5 are the results obtained in Okla- 
homa during April through July 1970. The measured 
velocities of the disturbances are compared with 
radar summaries and radar photographs. Spectra of 
the recorded waveforms are also presented. We 
discuss disturbances that produce S-shaped records 
and have periods in the acoustic domain, but 
have speeds appropriate to gravity waves. These 
disturbances suggest the presence of neutral winds 
in the F region. The approximate source power 
required for the acoustic disturbances detected in 
the ionosphere is discussed in section 6 along with 
the effects of the geomagnetic field on the detecta- 
bility of the waves. A microbarograph was installed 
in Oklahoma to see whether waves with periods such 
as those detected in the ionosphere were present 
near the surface. The (negative) results of the latter 
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experiment are discussed in section 7. Some statis- 
tical data on the occurrence of ionospheric distur- 
bances, the characteristics of the spectra, the influ- 
ence of geomagnetic disturbances, solar activity, 
etc., are presented in section 8. Finally, in section d 
we consider some of the possible causes of the dis- 
turbances, such as convective turbulence, shock 
waves from lightning, ionospheric absorption of 
sferics, and acoustic and buoyancy oscillations. 


2. HISTORY 


The idea that tropospheric weather affects the 
ionosphere began in the 1920s when Wilson (1925) 
was one of the first to suggest that the electric 
field of a thundercloud may cause ionization in the 
ionosphere. Other workers noted that reception of 
ionospheric radio echoes appeared to correlate with 
weather conditions. The appendix is a partial bibli- 
ography of publications in this field from 1925 
through 1971. 

Evidence for associations between weather and 
the ionosphere came from many parts of the world 
including the United States, the United Kingdom, 
Italy, and Australia. Regarding thunderstorms, it is 
interesting to recall the words of Watson-Watt 
(1933) in a discussion at the Royal Society of 
London. 

“I think it is clearly established by this series of 
observations that the sudden appearance of these 
bursts of abnormal ionization (in the E region) just 
precedes the arrival of thunderstorms in the area 
within 50 km of the recording station." 

Some of the correlations were misleading as, for 
example, the observation by Appleton et al. (1936) 
that *. . . abnormal E layer reflections occur most 
frequently in summer when local thunderstorms are 
at a maximum. . . ." We now know that this 
association results from two independent seasonal 
patterns. Associations between sporadic E increases 
and thunderstorms have also been found in India 
(Kesava Murthy 1963). 

Extensive studies were conducted in Australia by 
Martyn (1934) and his collaborators who found that 
day-to-day changes in the F region were associated 
with meteorological changes at the ground. Healey's 


(1936) calculations showed that lightning radiation 


could not account for the magnitude of the iono- 
spheric changes. 


After World War II, there were renewed attempts 
to establish correlations between ionospheric varia- 
tions and tropospheric weather, but with mixed 
success. Beynon and Brown (1951) inferred that the 
effects in both regions of the atmosphere were 
probably caused by solar activity. Evidence has been 
presented (e.g., Bauer 1958) for ionospheric effects 
of hurricanes. 

One of the main difficulties encountered in this 
work has been the inability of independent workers 
to verify the results. Indeed, in several cases, on 
the basis of more complete data, scientists have 
withdrawn earlier claims. 

That disturbances with acoustic periods in the 
lower atmosphere can couple into the ionosphere 
has been confirmed by observations of manmade 
explosions (Barry et al. 1966; Baker and Davies 1968) 
and earthquakes (Davies and Baker 1965; Yuen et al. 
1969). In November 1965, ionospheric wavelike 
disturbances with periods near 3 min were observed 
for several days over Boulder, Colo., and were 
thought to be associated with tropospheric phenom- 
ena (Davies 1966). Georges (1967) suggested that 
the relatively narrow spectrum of these disturbances 
was the result of the atmosphere filtering acoustic- 
gravity waves. Georges (1969) observed these 
3-min disturbances in the F region over the south- 
central part of the United States and made a major 
contribution by showing that they were associated 
with severe weather fronts passing through the 
observation area. Further studies of this phenomenon 
have been reported by Georges (1968) and by Baker 
and Davies (1969). They found that the following 
criteria must be fulfilled: 

1. Cloud tops must be greater than about 40,000 ft 
(12 km). 

2. Thunderstorm cells must lie within about 
200—250 km from the point on the ground directly 
below the ionospheric reflection point. 


Using the radar summary charts of the U.S. 
Weather Bureau (now the National Weather Serv- 
ice), Baker and Davies devised a severe weather 


‘index that ranged from 0 (no severe weather) to 5. 


They found that this index was always high (3 to 5) 
on days when disturbances were observed in the F2 
region. However, on days of high severe weather 
index, ionospheric disturbances were not always 
observed. Similar observations of ionospheric dis- 
turbances associated with thunderstorms in Ala- 


bama have been made by Detert (1969). 
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3. IONOSPHERIC RADIO SOUNDING 
A. Magneto-Ionic Theory 


To understand the radio sounding techniques 
used in this investigation, we shall briefly review 
the propagation of radio waves in a cold plasma con- 
taining a magnetic field. More detailed theory is 
given elsewhere (e.g., Ratcliffe 1959; Budden 1961; 
and Davies 1969). The interested reader should 
consult these books. 

When a radio wave enters the ionosphere, it is 
split into two oppositely polarized waves. On radio 
frequencies above the electron cyclotron or gyro- 
frequency (fu ~ 1.4 MHz), the left circular polarized 
wave (looking along the direction of the earth's 
magnetic field) is called the ordinary wave. The right 
circular polarized wave is called the extraordinary 
wave. On a single frequency, there are, in general, 
two radio echoes reflected from different heights. 
Because of collisions between electrons and neu- 
tral molecules, energy is absorbed from the wave. 
Most of the absorption occurs in the D region (60 to 
90 km) where the collision frequency v is high. The 
absorption of the extraordinary wave is usually 
greater than that of the ordinary wave; thus, the 
latter wave usually dominates the composite signal. 

The maximum frequency reflected from a par- 
ticular layer of the ionosphere is called the critical 
frequency of the layer. The ordinary wave critical 
frequency f, is related to the maximum electron 
density Nm of the layer by the expression 


Nm= 1.24 X 1019/2 m-3 (1) 


where f, is in megahertz. The critical frequency of 
the extraordinary wave f; is related to f, by 


fo Ce CH (2) 
Hence, f; > fo. 
The complex refractive index of the ordinary (+) 
and extraordinary (—) waves of frequency f are 
given by the Appleton equation 


M? — (y — ix)? -1—XJD (3) 


where 


D—1-—iZ —[Y32(1 — X — iZ)] 


+ V[IY?94(—X-iZ) FY? 


474-478 O - 73 - 2 


and 
X = (frl f)? 
where 
fn= plasma frequency, 
Y —fylf, Y, —Y cos 0, Y,—Y sin 6, 
Z —v[2mf, 
and 
i= v—]. 


Here, 0 is the angle between the direction of wave 
propagation and the geomagnetic field. 

The important feature to notice in eq (3) is that 
the propagation depends on wave frequency 
(dispersion) and on direction of propagation (anisot- 
ropy) In dispersive propagation, the wave packet 
(energy) travels with the group speed that is, in 
general, different from the phase speed (v) with 
which the individual wave components travel. Hence, 
the individual waves tend to pass through the packet. 
In anisotropic propagation, the direction of motion of 
the packet (ray direction) differs from that of the 
phase normal The component u of the group 
velocity V, in the direction of the wave normal (k) 
is given by 


Qo do) 
pe Vgsec ced uen (4) 


where o is the angular wave frequency and c is the 
free space speed; o, the angle between the wave 
normal and the ray, is given by 


` Län 
Ge (5) 


Ít is important to bear these properties in mind 
because, as we shall see later, they apply also to the 
propagation of atmospheric waves. 

Because of the anisotropy of the ionosphere, for a 
radio wave incident vertically on the ionosphere (in 
the Northern Hemisphere), the ordinary wave is 
reflected to the north of the point of entry; and the 
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FIGURE 1. — Sample Doppler record (9.9 MHz, Havana, Ill., to Boulder, Colo., 13-14 June 1967) 


extraordinary wave is reflected to the south. Hence, 
a disturbance traveling in the ionosphere will, in 
general, affect one wave before the other (Davies 
1969, ch. 7). 4 

The sweep frequency radar or ionosonde produces 
records of time delay (or virtual height) versus 
frequency for both ordinary and extraordinary 
waves. From these records, it is relatively easy to 
obtain the electron density profile up to the peak of 
the F2 layer (Wright and Smith 1967). 


B. Doppler Technique 


When a change takes place in an electron density 
profile, the total phase path P of a radio signal 
changes; and the instantaneous frequency shift 
Af of a radio echo of frequency f is given by 


A f dP 
Kcd (6) 
where 

P= f. u cos a ds. (7) 


Here, ds is an element of the ray path s. As seen 
from eq (3) and (5), u and o are rather complicated 
functions of frequency, direction, electron density, 
etc. In broad terms, however, two types of change 
predominate: 

1. A change of path length s. 

2. À change of refractive index at levels below 
reflection. ( 
These two types of phase path change (provided we 
neglect collisions and the geomagnetic field) can 
be distinguished by observations on several close 
frequencies. In type 1, 


Af f. (8) 
In type 2, 
A f c: 1/f. (9) 


For the disturbances discussed in this paper, the 
Doppler shifts obey eq (8); henceforth, we shall 
assume that the observed shifts result from changing 
heights of reflection. 

The continuous measurement of the Doppler 
shift is carried out by mixing the ionospheric echo 
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with a local offset reference signal and recording 
the beat frequency on magnetic tape moving 
l in./min. The frequency analysis is accomplished 
by fast playback of the tape, resulting in a frequency 
multiplication of about 1800. The audio signal thus 
produced is analyzed by conventional means (Davies 
1969, sec. 9.3) and displayed on a facsimile chart. 
A sample record is shown in figure 1; from this, we 
see that the Doppler shifts are about 1 Hz. The 
smallest frequency shift that can be measured 
with this technique is about 0.1 Hz. 


C. Circuits 


The radio paths used in this paper are shown in 
figure 2, together with the sounding frequencies 
used and the path midpoints. We dwell in detail on 
the spaced transmitter experiment in Oklahoma 
during the summer of 1970; we use the remaining 
data in an auxiliary role. Oklahoma was chosen for 
the observations because this State has the highest 
incidence of severe storms, especially tornadoes, 
per unit area of all States (National Weather Service 
1970). Furthermore, weather radar data were avail- 
able. The particular circuits shown in figure 2B 
were chosen partly because facilities were made 
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available by NOAA's National Severe Storms Lab- 
oratory at Norman and by the Electrical Engineering 
Department of Oklahoma State University at Still- 
water. A separation of about 50 km between iono- 
spheric reflection points at the corners of an 
approximately equilateral triangle was selected 
because of the success of a somewhat similar experi- 
ment near Boulder (Davies and Jones 1971a). 
Ultimately, the distances selected were small 
because the horizontal trace velocities involved 
were about 1 km/s and the resulting time displace- 
ments less than a minute. Of course, if the separa- 
tions are increased, the cross correlation between 
the various displaced Doppler curves is decreased. 
The choice of radio frequencies was restricted by 
ionospheric conditions to the range of frequencies 
that reaches, but does not penetrate, the F2 layer. 


D. Data Analysis 


With the arrangement of spaced transmitters on 
3.3, 4.0, and 5.1 MHz used in Oklahoma during 
April to July 1970, nine simultaneous Af (t) records 
were obtained. The beat frequencies were set about 
1 Hz apart to distinguish the three traces on any 
one carrier. The velocity can be determined from 


STILLWATER, OKLA. 
3.34, 4.05 5.1- MHz 
Transmitters 


EL RENO, OKLA. 
3.34 4.07 5.1- MHz 
Transmitters 


SHAWNEE, OKLA. 
3.3, 4.04 5.1- MHz 
Transmitters 


NORMAN, OKLA. 
3.34, 4.04, 5.1- MHz 
Receivers 


FIGURE 2. — (A) radar summary chart of the radio paths from Longbranch to Boulder and Norman to Stillwater; (B) spaced transmitters 
in Oklahoma 
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FIGURE 3. — Horizontal trace velocity Vz from time displacements 


the time displacements by a procedure similar to 
that described by Davies and Jones (1971a), on the 
assumption that the disturbance is plane. This 
assumption is examined in subsection 4E. 

The horizontal trace velocity Nu is obtained by 
the construction shown in figure 3 that represents 
a plane wave front moving across the three (reflec- 
tion) points E, U, and S. Let Ter and Tee be the 
time delays between the patterns at E and U and 
at E and S, respectively; then the apparent or 
trace velocities Nor and Vgs are given by 


Vev= 56.0 km/Tzgvu 
and (10) 


The magnitude and direction of Nu can be found 
from the following equations. The direction o (east 
of north) is given by 


Vev sin Ugu — Ves sin bes | 


: 11 
Vru COS Urku— Ves COS Wes ( ) 


cot g=- 


the magnitude Vy is given by 


FIGURE 4.— Derivation of true velocity V from horizontal and 
vertical trace velocities No and Vy 


sin? (Wes— ev) . 1 , l1 2 cos(Wzs— ev) 


d Bu Vis 
(12) 


Here, Wev and Wes are the angles between geo- 
graphic north and the lines EU and ES, respectively. 

To obtain the true velocity V, we also need the 
vertical trace velocity Vy found from the time 
displacements of the radio signals on different 
frequencies and the heights of reflection of the 
radio signals. The reflection points for a given 
transmitter and receiver location are assumed to 
lie on a vertical line. Errors introduced by this 
assumption are discussed in subsection 4E. The 
heights of reflection can be found by true height 
analysis of ionograms (Wright and Smith 1967). 
Knowing the trace velocities Nu and Vy, we 
obtained the true velocity V from the construction 
in figure 4 showing that the true speed is smaller 
than the vertical and horizontal trace speeds, 
which are, therefore, not components of the true 
velocity. 

The time displacements can be determined in 
two ways. One method is to locate similar features 
on the records (e.g., peaks, troughs, etc.) and to 
average the time displacements over as large a 
sample of cycles as possible. In the other method, 
the traces are sampled at intervals, short compared 
with a cycle period; and the time displacements are 
measured from the cross-correlation functions 
(Jones 1969). 

The application of these two methods to the 
disturbance of 0425-0505 UT on June 11, 1970, is 
illustrated in figures 5 and 6. From the positions of 
the peaks and troughs, the average time displace- 
ments are 7£5—35.8 s and Tey=— 23.3 s. The 
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Shawnee, Okla. (Upper Trace) 
EI Reno, Okla. (Middle Trace) t to Norman, Okla. 5.1 MHz 
Stillwater, Okla. (Lower Trace) 


VARIA 


gr 0430 0435, wa ou 0490 0455 5 05 0505 


FIGURE 5.—Method of digitizing the Doppler records for the 
disturbance of June 11, 1970, at 0425—0505 UT 


Doppler record has been carefully scaled every 30 s 
(fig. 5) to find the time displacement by cross correla- 
tion. The instantaneous Doppler frequency is repre- 
sented by the distance d on the record from the 
traces to the reference line. The reference line 
follows the long-period trends of the traces to 
eliminate their influence on the cross correlograms. 
The cross-correlation curves (fig. 6) give Tgs = 38 s, 
Tgu— —22 s, and Tys=58 s. Similar analyses 
have been carried out for the disturbances of April 
30 and May 13, 1970. The magnitude Ho and azimuth 
o obtained by these two methods are given in table 1. 

In general, the agreement is good except for the 
azimuth on April 30 when the trace velocity is 
large, which implies that the direction of propaga- 
tion is nearly vertical. 


4. ACCOUSTIC-GRAVITY WAVES 


A. Refractive Index and Parcel Orbits 


Acoustic-gravity waves are oscillatory motions 
traveling in a stratified atmosphere in which the 


R 


0425-0505 UT 


Cross - Correlation Coefficient 


=3 E: " 0 1 2 3 


FIGURE 6.—Cross-correlation functions for the records from 
El Reno, Shawnee, and Stillwater on June 11, 1970 


TABLE l.— Comparison of horizontal velocities 
obtained by similar features and cross correlation 


Date Vu c 
(1970) UT (m/s)  (deg.) 
Peak and trough 850 158 
am — ee, correlation 840 157 
Peak and trough 2100 54 
3 —1040 
co — we correlation 2600 18 
Peak and trough 900 172 
—0 
i —-— rs correlation 880 176 


forces of buoyancy are comparable to the compres- 
sional forces. Detailed studies of the propagation of 
these waves in nonuniform absorbing and windy 
atmospheres have been given by a number of 
authors (e.g., Hines 1960; Midgley and Liemohn 


1966; and Pitteway and Hines 1963 and 1965). 


The propagation of acoustic-gravity waves in a 
temperature and density stratified and flat atmos- 
phere is dominated by certain characteristic 
frequencies. Neglecting the earth's rotation, there 
are two such frequencies: the acoustic cutoff fre- 
quency «4 and the buoyancy (also called the 
Brunt-Vaisala) frequency op given by 


WA = Wa (13) 
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Acoustic 
Waves 


Gravity 
Waves 


FORBIDDEN — 


Wg Wo 


FIGURE 7.—Frequency ranges of acoustic and gravity waves 


and 
ld 
Y x 
of = (173 m) (14) 


where y is the ratio of the specific heats at constant 
pressure and constant volume and H' is the gradient 
of scale height H. The frequencies wa and wg are 
the isothermal values of the acoustic cutoff and 
buoyancy frequencies given by 


wa = p (15) 
and 
oj pl. (16 
Here, 
C? = yg H (17) 


is the square of the sound speed, and g is the 
gravitational acceleration. 

The physical significance of these two frequencies 
is as follows. If we squeeze the whole atmosphere 
vertically and then release it, the atmosphere will 
oscillate vertically about its equilibrium position at 
wa. On the other hand, if we displace a parcel of 
the atmosphere vertically and then release it, the 
parcel will oscillate at wg. The two frequencies 


are related by 


fwb ` a. (18) 


2 
@ 

a 2 
y 9 


Hence for an atmosphere with diatomic molecules 
for which y = 1.4, 
Wg X 0.9 Qa. (19) 


‘From Davies et al. (1969), the refractive index u of 


an acoustic-gravity wave of angular frequency w 
propagating at an angle ¢ to the vertical in an iso- 
thermal atmosphere is given by 


Been DONT SU DENN 
CUPS p ] — Y? sin? $ 


(20) 


where 


V is the phase velocity, and Eis the wave propaga- 
tion parameter. 

Equation (20) reveals the frequency dependence 
(dispersion) and the angular dependence (anisot- 
ropy) of the propagation. It also shows that there 
are two families of waves that can propagate. One 
family has X (and, hence, Y) less than unity. These 
high-frequency waves constitute the acoustic 
domain. The low-frequency waves have X and 
Y sin @ greater than unity and are called gravity 
waves. Thus, acoustic waves exist at frequencies 
greater than the acoustic cutoff frequency, whereas 
gravity waves exist at frequencies less than the 
buoyancy frequency. In between, there is a for- 
bidden frequency range as sketched in figure 7. 
The u(ġ) curves, parametric in Y, are shown in 
figure 8 for acoustic waves and figure 9 for gravity 
waves, together with the particle orbits given by 


u sin $ 


PET TPA] {—u cos $ + i (X — Y2)12). 


(21) 


E 
d 


Here, £ and ¢ are the horizontal and vertical displace- 
ments (or velocities), respectively; and i= V—1. 
The orbits are ellipses that, under certain con- 
ditions, become circular or linear. For acoustic 
waves, the orbit is circular for horizontal propaga- 
tion at a frequency wa V2/y. The orbits are linear 
with vertical propagation and almost the same at 
frequencies approaching the acoustic cutoff fre- 
quency for all directions of propagation. As the wave 
frequency departs farther from and above the 
acoustic cutoff frequency, the orbits become 
essentially longitudinal. 
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R=0.09 
y=90° 


R=0735 R=| 
y=90° 


FIGURE 8. — Refractive index curves (u surfaces) and parcel orbits for acoustic waves in an isothermal atmosphere (y= 1.4). R is the 
axial ratio of the orbit, and y is the inclination of the major axis from the horizontal. 


For gravity waves, propagation is limited to 
certain directions. The minimum and maximum 
values of $ are 


i o 
min = sin^! (>) 


Wg 


and (22) 
max = 7 — sin! (>) 


Wy 


The sense of rotation of the air parcels for gravity 
waves is the reverse of that for acoustic waves. 
The orbit is circular for horizontal propagation on 
a frequency w=Vy/2 wg. As $ approaches its 
asymptotic values, the parcel orbits become linear 
and transverse to the direction of propagation. 
Some of the properties of acoustic waves and 


gravity waves are compared and contrasted in 
table 2. 
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Y =1.05 


FIGURE 9. — Refractive index curves (u surfaces) and parcel orbits for gravity waves in an isothermal atmosphere (y= 1.4). R is the axial 


ratio of the orbit, and y is the inclination of the major axis from the horizontal. 


TABLE 2.—Some properties of acoustic and gravity B. Ray Tracing 


waves in an isothermal atmosphere 


We have seen that the propagation of acoustic- 
Parameter Acoustic waves Gravity waves gravity waves is dispersive and anisotropic. As 
mentioned in subsection 3A, the flow of energy 
differs in both direction and speed from the flow 
Y= olo (Y «1 (Ysinó 21 of phase. The component u of group velocity V, in 


X = (val o)? (X«1 (X»1 


Refractive index {p<l WIEN the direction of wave propagation is obtained by 
Phase speed {v>C Keep substituting eq (20) into eq (4). This gives 
— &in-1 > 
Propagation angle {0< ġ<2r E - moe V(w?— w2) 
2 sin^! (1/Y ) u= KN COS a= ———— t : 
H 2 L3. t 
Air parcelrotation ` ( Clockwise { Anticlockwise o bo, sin dh 


Energy flows up when 
phase travels down 
and vice versa. 


| Same quadrant 


Energy vector S : 
as propagation 


The angle o between the wave normal and the ray is 


peas Cleat ee 
1 —p?Y? sin? $ 


11 
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Z Ü 


FIGURE 11.— Relative directions of group propagation without 
wind (V) and with wind (S) 


likewise obtained from eq (2) and (5): 


. Y? sing cos o. 

NUT ] — Y? sin? $ 
Equation (24) shows that, for acoustic waves, 
a (20) lies in the same quadrant as ¢; but for 
gravity waves, a(« 0) is in a different quadrant 
from $. This means that, for upward progression 
of energy flow in a gravity wave, the phase normal 
points downward. Applying eq (23) and (24) to the 
propagation of gravity waves illustrates the beaming 
properties of the propagation. Figure 10 (p. 11) shows 
surfaces reached in unit time by gravity-wave 
packets radiated from a (hypothetical) point source. 
The limiting angles of the energy propagation are 
given by 


tan (h — Q) mar — VY?—1. 


min 


(25) 


Let us now derive the equations necessary to 
trace acoustic-gravity waves in a slowly varying 
atmosphere in the presence of a horizontal wind U in 
the plane of propagation. 

The refraction condition for these waves in a 
windy atmosphere is expressed by the constancy of 
the horizontal trace velocity Vu: 


m.s +U=Vu (26) 
sin dh 

(see, e.g., Lamb 1960, sec. 75). We start at some 
reference level in the atmosphere with an initial 
value of $. To find Vy (for a known U), one must 
know V. This can, in principle, be found from eq 
(20); but the latter applies only to an observer at 
rest with respect to the medium. If w is the wave 
frequency observed on the ground, the intrinsic 


(24) 


frequency Q as seen by the observer moving with 
the wind is 


Q—o(1-— UINal (27) 


(Hines and Reddy 1967); €) must be used in eq (20) 
in place of w. Combining eq (20), (26), and (27) 


yields a quartic that can be solved for Vu. The ray 


direction B is obtained directly from the geometry 
of figure 11: 


EE? V, cos(ġ — a) (28) 
From eg (23), 
tan ($ — o) = (1 — w?/”) tan dh, (29) 


Hence, the horizontal coordinate x of a point on the 
ray at a height z is given by 


x= | "ran B-dz. (30) 


20 


The time of flight t is given by 


" dz 
pel eee’ 31 
J o Va cos(ó — a) Ge 
The phase, relative to that at the reference height 
Zo, is given by 


o=| CP. cos asec B dz. (32) 
Zu 


C. Ray Paths in Model Atmosphere 


The equations developed above have been inte- 
grated by a computer for the 1962 U.S. Standard 
Atmosphere (COESA 1962). The height variations 
of molecular temperature Ty, sound speed C, and 
the characteristic periods Tg, ts, and Ta for this 
atmosphere are given in figures 12, 13, and 14. 

Because the speed of sound increases above about 
90 km, acoustic-gravity waves are refracted (and may ` 
be reflected). For acoustic waves, there exists an 
iris in this region of the atmosphere through which 
waves can penetrate to the F region. Gravity 
waves, which cannot propagate at angles to the 
vertical less than the critical angle c, can only 
reach the F region through an annulus. This latter 
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FIGURE 12.— Height variation of molecular temperature Ty 
for the 1962 U.S. Standard Atmosphere 
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FIGURE 13.— Height variation of sound speed C for the 1962 
U.S. Standard Atmosphere 


property results in an angular filtering of gravity 
waves propagating into the F region from sources in 
the lower atmosphere. Hence, F-region traveling 
disturbances originating in the lower atmosphere 
will propagate only at angles near the asymptotic 
value. This is illustrated in figure 15. Recall that, 


Height, km 


Period, min 


FIGURE 14.— Height variations of acoustic cutoff and buoyancy 
periods for the 1962 U.S. Standard Atmosphere 


I962 U.S. STANDARD ATMOSPHERE 
T 730 min 
Source Height = IO km 


90 120 150 180 
Propagation Angle 4, aeg. 


FIGURE 15.— Height variations of the propagation angles of 
30-min gravity waves from a point source at 10 km (from 


Chang 1969) 
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FIGURE 16.— Acoustic ray paths in the 1962 U.S. Standard Atmosphere from a source at 15 km and with periods of 3.5 and 4.0 min 
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FIGURE 17. — Variation of ground range with the initial propaga- 
tion angle 
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FIGURE 18.— Variation of reflection height with the initial 
propagation angle 
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1962 U.S. Standard Atmosphere 
T=O.5min, Source at I3km 


1962 U.S. Standard Atmosphere 
T= 3.0min, Source at 20km 
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FIGURE 19.— Acoustic ray paths (from an elevated source) showing ducts 
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FIGURE 20. — Variation of ground range with the propagation angle $, measured at the level of radio reflection, for a source at a 15-km 
height and with periods of 3.5 and 4.0 min 
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for upward energy propagation, the propagation 
vector is directed downward. 

The equations in subsection 4B have been 
integrated to give ray paths of acoustic waves with 
periods of 3.5 min and 4.0 min for no winds and for a 
source height of 15 km. These paths (fig. 16, p. 15) 
show that, above 100 km, the iris changes with wave 
period. Figure 16 shows that, for a given height in the 
atmosphere (e.g., 100 km), the aperture through 
which the waves flow for a given initial angle ġo de- 
creases with increase of wave period. However, 
when the wave period exceeds the minimum value of 
the acoustic cutoff period (which usually occurs near 
85 km), total internal reflection occurs at all angles. 
Under these conditions, energy can only “tunnel” 
through the barrier, and "full wave" theory must 
replace ray tracing. It is interesting to note that, 
near the 85-km level, the propagation angle œ for a 
given initial angle ġo tends to increase because the 
period is approaching the acoustic cutoff period; 
but this may be more than offset by the decrease 
brought about by the slow sound speed. 


Notice also that the initial propagation angle 
required for reflection at a given height (around 
150 km) increases as the period increases. Figure 17 
shows how the ground range varies with initial 
propagation angle for waves of different periods. 
For waves with periods less than about 4.0 min, there 
are two initial angles for a given range greater than a 
minimum range. This property is reminiscent of the 
oblique propagation of ionospheric radio waves 
(Davies 1969, ch. 12) in which the minimum distance 
is called the skip distance and the rays with the 
smaller and larger values of ¢ are called high-angle 
and low-angle rays, respectively (high refers to the 
elevation angle). For waves with periods between 4 
and 5 min, there are discontinuities at the angles at 
which the wave penetrates a layer formed by the 
temperature variations shown in figure 12. The 
heights of reflection also change abruptly at these 
angles, as can be seen in figure 18. The reflected 
waves are trapped between the reflection levels (one 
of which may be at the ground) and propagate as in a 
waveguide (fig. 19). 

As we shall see later, our main concern in this 
paper is with acoustic waves having periods in the 
approximate range of 2 to 10 min; these waves pene- 
trate to the F region. Hence, we shall no longer 
concern ourselves with reflected waves. 


In practice, we measure the magnitude and direc- 
tion of the phase speed in the F2 region at the radio 
reflection levels. From a knowledge of the propaga- 
tion angle and the height of reflection (from iono- 
gram analysis), the ground range to the source can 
be calculated as a function of $ for various reflection 
heights. Families of curves are shown in figure 20 
for waves with periods of 3.5 and 4.0 min giving the 
variation of ground range with angle in the ionos- 
phere for various heights. Comparison of the two 
families of curves in figure 20 shows very little 
difference. This is important because, as we shall 
see, disturbances associated with thunderstorms 
contain a range of periods. The reason for the close 
similarities in the curves for the two periods is that 
most of the ray refraction occurs above 100 km 
where the acoustic cutoff and buoyancy periods are 
considerably greater than the wave periods (i.e., 
X<<1, Y« «]1) therefore, the propagation is 
essentially nondispersive. Hence, the major error in 
these ground range calculations is caused by un- 


certainties in the temperature model above 100 km. 
Figure 21 shows curves of constant phase for 


waves of period 3.5 min emanating from a point 
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FIGURE 21.— Surfaces of constant phase in the 1962 U.S. 
Standard Atmosphere (no wind) for waves of period 3.5 min 
emanating from a source at 15 km 
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source 15 km above ground. The interesting feature 
to note here is that, for the heights of concern to us, 
there are only about two wavelengths between the 
source and the area of observation in the ionosphere. 
Under these conditions, it is clear that ray theory 
cannot be used to accurately determine the wave 
properties because the wavelengths are much larger 
than the scale heights of the medium. 


D. Determination of Velocity of 
Ionospheric Disturbances 


The determination of the three-dimensional 
velocity of the ionospheric disturbance requires the 
horizontal and vertical trace velocities. The hori- 
zontal trace velocity is determined by the method 


described in subsection 3D. To find the vertical trace . 


velocity, we need to know the locations of the 
reflection points of the radio waves of different fre- 
quencies and the appropriate time displacements. 
Unfortunately, in the present program, height infor- 
mation at the time of the measurements was not 
always available; in several cases, the amplitudes 
of the disturbances on the lower frequencies of 
3.3 and 4.0 MHz were too small to permit measure- 
ment of the time displacements. However, we are 
satisfied that, for all waves with acoustic periods 
discussed in this paper, the phase propagation was 
upward. 
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FIGURE 22.— Variation of propagation angle œ with horizontal 
trace velocity Nu 


To obtain the propagation angle @ at any height 
with no observations of the vertical trace velocity, 
we proceed as follows. From the model atmosphere 
given in figure 13, we have the sound speed C and 
the characteristic periods Tq and Tg. In the absence 
of winds, eq (26) gives 


sin d = $ (33) 


Substituting eq (33) into eq (20), we obtain 


4. i-um T 
y C? riy 


(34) 


Thus, for a given atmosphere, V and Vy are unique 
functions of $. The variation of $ with Vy is illus- 
trated in figure 22. 


E. Some Possible Errors in the Velocities 


It is now appropriate to discuss some errors that 
the assumptions introduce into the velocities. One 
of the important assumptions is that the wave front 
is plane, which will not be true for acoustic waves 
from severe thunderstorms. From observation, we 
know that the sources are about 300 to 500 km from 
the center of the triangle of sounding points. To get 
an idea of the error introduced by the curvature of 
the wave front, we assume a radius of curvature of 
400 km and a 50-km equilateral triangle. The maxi- 
mum error in direction occurs when propagation is 
parallel to a side of the triangle and amounts to 
0.96°. There is no error in direction when propaga- 
tion is parallel to a median. By contrast, the error 
in the magnitude of the velocity is a maximum for 
propagation along a median, 1.8% in this case. 

Another assumption is that the reflection points 
for 3.3, 4.0, and 5.1 MHz are stacked vertically. 
The fact that they are not (because of magneto-ionic 
deviation) causes errors in the determination of the 
vertical trace velocity. Because the deviations, 
which lie in the magnetic meridians, differ for 
different frequencies, the triangles formed by the 
reflection points will be displaced and distorted. 
The midpoint triangles are shown in figure 23 for 
the concentric density profile shown in figure 24 
(top). Electron density profiles were unavailable in 
Oklahoma, so we have used a Boulder, Colo., pro- 
file to illustrate these errors. 
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OKLAHOMA DOPPLER, REFLECTION POINT GEOMETRY 
N(h) PROFILE BOULDER, APR. 25, 1968, 0100 UT 
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———— Actual Triangle at 5.1 MHz 
—— —— Actual Triangle at 4.0MHz 
HERE Actual Triangle ot 3.3 MHz 
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FIGURE 23.—Triangles formed by midpoints of radio paths and 
by the reflection points on 3.3, 4.0, and 5.0 MHz 


The neglect of the lateral deviations of the 
ordinary waves shown in figure 24 (bottom) can lead 
to appreciable errors in the vertical trace speed. 
This can easily be seen qualitatively by considering 
a vertical wave front moving south. It will affect the 
3.3, 4.0, and 5.1 MHz successively. On the assump- 
tion of vertical stacking of the reflection points, it 
would be deduced that the wave has an upward 
component of velocity. The error in tilt of the 
propagation vector with respect to the vertical is 
shown in figure 25 for various azimuths o. The 
maximum error of 9.3° occurs for meridional 
propagation and decreases to zero with propagation 
from west to east and east to west. The error in the 
true speed is shown in figure 26; note that the 
maximum errors of about 8.5% occur for propaga- 
tion angles near 45° and near 135° with the vertical. 
However, as pointed out by Davies and Jones 
(1971a), the most important source of error is a 
mistake in distinguishing between an ordinary 
and an extraordinary trace. 
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FIGURE 24. — Model I (M-I) ionogram (top) and the lateral devia- 
tions of the reflection points of the ordinary waves (bottom) 
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FIGURE 25. —Errors in speed (as function of propagation angle dh 
and azimuth ol caused by lateral deviation 
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5. OBSERVATIONS 


ie d Hel A. General Remarks 
OIOO UT 
In this section, we present the experimental 
evidence, gathered during the spring and summer 
of 1970, for wavelike disturbances in the ionosphere 
and their relationship to thunderstorms. In par- 
ticular, we use the theory developed in the preced- 
ing sections to try to locate the sources of the 
disturbances. The source locations are based on a 


source height of 15 km. 


B. April 18, 1970 


The Doppler record for this event is shown for a 
radio frequency of 5.1 MHz (fig. 27). The Doppler 


shifts on the lower frequencies were too small to 


V' = Apporent Speed Neglecting M-I Deviation 


SE be usable. This arises from the higher reflection 
we 30 60 90 120 150 m height of the 5.1-MHz signal and the resulting 
Send greater amplification of the acoustic wave because 
of the lower gas density. 
FIGURE 26.— Errors in the propagation angle caused by lateral The Doppler traces in figure 27 have been sep- 


deviati a function of the angle) i 
— ER arated by setting the three transmitters at fre- 


quencies about 1 Hz apart. The displaced lines of 
dots near the bottom of the record provide a time 
code. 


SHAWNEE, OKLA. (UPPER TRACE) 
EL RENO, OKLA. (MIDDLE TRACE TO NORMAN, OKLA. 5.1 MHz 


STILLWATER, OKLA. (LOWER TRACE) 


E. LE 24.5... | en d 


z 
a 
N 


0700 0800 0900 


FIGURE 27.— Doppler record (5.1 MHz) for the disturbance of Apr. 18, 1970, at 0500—0900 UT 
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The locations and heights of the tops of the 
thunderstorm cells, as obtained from the National 
Weather Service's radar summary charts, are 
shown in figure 28. The horizontal trace velocities 
are shown by arrows, and the points of origin deter- 
mined by ray tracing are shown by the symbol C. 
Extrapolating in the direction opposite to Vy (fig. 
28A), we come to the cells designated 480 and 510 
(hundreds of ft) These thunder cells meet the 
criteria (established by Baker and Davies 1969 and 
by Davies and Jones 19715) that the storm must lie 
within a radius of about 200 km of the sounder and 
have cloud tops higher than 40,000 ft. We conclude 
that one or both of these thunder cells is responsible 
for the acoustic wave effects detected in the 
ionosphere. 

One of the more interesting features of this storm 
is that, from our ionospheric soundings, we could 
follow its movement to the northeast. This can be 
seen from the sequence of radar summaries in 
figure 28 and on the radar photographs of figure 29. 
We see that there is good overall agreement 
between the position of the storm and the direction 
of the trace velocity. The numerical data for this 
and the following disturbances are shown in table 3. 
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TABLE 3.— Summary of data from acoustic wave 


events of 1970 

Vu c 
Date UT (m/s) (deg.) 
Api 1B uus essa ean nto MOUSE 2750 97 
|i o ` E TEEN am os sc 1500 126 
AME EE ORUM 1800 147 
E ENT, ET, 850 206 
An MAL ecol EE MEL 2700 54 
AT We OSCR ERG 900 172 
June IIS EN VE er 850 158 


SHAWNEE, OKLA. (UPPER TRACE) 
EL RENO, OKLA. (MIDDLE TRACE 
STILLWATER, OKLA. (LOWER TRACE) 


TABLE 4. — True reflection heights (at the Oklahoma 
experimental site) estimated from White Sands, 
N. Mex., ionograms 


Date True heights (km) 
(1970) UT 
3.3 MHz 4.0 MHz 5.1 MHz 

Asp MEL S une 0000 299 312 352 
Age 35 uelle. 2800 106 130 174 
Apo dE. aes... 1900 112 131 183 
Apt JN Lun oos 2000 134 162 207 
Ar. 20,0951 7- M00 333 352  Áf,«5.1 MHz 
Meggie, 70000 144 174 212 
May.3............... "0200 230 236 249 
May 153. 5ue.e. o 9900 306 326 313 
June 1l. lassus OE 256 266 284 


TO NORMAN, OKLA. 5.1 MHz 
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FIGURE 30.—Doppler record for Apr. 29, 1970, 2100 UT to Apr. 30, 1970, 1200 UT 
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FIGURE 31.— Oklahoma City weather radar photographs (250-n.mi. range) and horizontal trace velocities (V; ) for Apr. 30, 1970 
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FiGURE 32. — Radar summaries with horizontal trace velocities for Apr. 30, 1970 


Unfortunately, ionograms were not available from City. The results are summarized in table 4 (p. 23). 
Oklahoma. In lieu of these, the true heights of The times listed here are 0.5 hr earlier than those 
reflection (needed to determine the ground ranges) when the ionograms were obtained, to allow for the 
were obtained by analysis of ionograms from White difference in local times. 


Sands, N. Mex., about 800 km west of Oklahoma 
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FIGURE 33. — Doppler record for May 12, 1970, 2300 UT to May 13, 1970, 0500 UT 


C. April 29-30, 1970 


This disturbance shown in figure 30 (p. 23) is of 
particular interest because, near 0830 UT, a tornado 
accompanying it passed through the triangle formed 
by our transmitters. During the afternoon of April 
29, a disturbance producing an S-shaped record 
occurred, which is discussed in subsection 5G. 
Evidence of shorter (3- to 5-min period) waves 
appears near 0030 UT but does not become promi- 
nent until about 0430 UT. At this time, unfortun- 
ately, the ionospheric signals are spread because of 
reflections from small ionospheric irregularities 
associated with enhanced geomagnetic activity. 
Hence, the record is unsuitable for detailed analysis 
until shortly before 1000 UT. Velocities, calculated 
from time differences near 1000 and near 1030 UT, 
are superimposed on the radar films in figure 31. 
There is an indication, which we believe to be real, 
that the source of the disturbance changes (rather 


abruptly) from one source within the storm system 
to another, as shown in the radar summaries of 
figure 32. Near 0945 UT, the primary source is 
probably that marked A in the left radar summary in 
figure 32. By 1045 UT, the primary source was 
either B or C (or possibly both) in the right sum- 
mary. All these cells have the properties necessary 
for producing observable ionospheric disturbances. 


D. May 12-13, 1970 


The Doppler record (fig. 33) shows an S-shaped 
disturbance followed by a 3- to 5-min wave. The 
splitting of the traces is caused by ordinary and 
extraordinary waves having different reflection 
points. The radar summary in figure 34 shows no 
activity over Oklahoma near 0245 UT, but there are 
two small cells in northern Texas associated with 
tornadoes and two large cells over northeastern 
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FIGURE 35. — Oklahoma City weather radar photograph (250-n.mi. 
range) on May 13, 1970 at 0359 UT. V,-—900 m/s, and 
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FIGURE 34.— Radar summary with horizontal trace velocity for 
May 13, 1970 
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FIGURE 36. — Doppler record for June 11, 1970, at 0300—0700 UT 


Kenneth Davies and John E. Jones 27 


100° W 90° W 


Vu = 850 m/s 
c 7158? 


} 0425 -O505 UT | 


40°N 


o 100 200 300 
EL 


km JUNEI!, 1970,0445UT 


FIGURE 37. — Radar summary for June 11, 1970 


Kansas. The horizontal trace velocity suggests 
that the source (or sources) of the ionospheric 
disturbance is in the Kansas storm. Of course, 
the waves may be coming from several sources; and 
this may well be the reason for the “ragged” appear- 
ance of the disturbance. The radar photograph in 
figure 35 shows the two centers of activity but is not 
nearly as impressive as the radar summaries. This 
event is one of the rare exceptions to the rule 
that the thunderstorm must lie within a radius of 
about 200-250 km of the sounding area to produce a 
detectable ionospheric disturbance. 


E. June 11,1970 


The Doppler record for this event is shown in 
figure 36. The radar data (figs. 37 and 38) show a 
single storm in the region; the direction of the 
velocity vector agrees well with the weather data. 


YR" 


FIGURE 38. — Oklahoma City weather radar photograph (250-n.mi. 


range) with the horizontal trace velocity for the disturbance 
of June 11, 1970, at 0434 UT. V, — 850 m/s 


F. Power Spectra of Acoustic Disturbances 


From a physical viewpoint, an important charac- 
teristic of these ionospheric disturbances is the 
power spectrum. The "power" in the Doppler record 
is proportional to the mean squared Doppler shift 
and is therefore porportional to the mean squared 
velocity of the air parcels excited by the acoustic 
wave. Thus, the power spectrum of the Doppler 
record is also a power spectrum of the wave itself. 
The records were scaled every minute, and the 
spectra were obtained by a program written by 
Lewis (1969) based on the method of Welch (1967). 
Before discussing the actual spectra obtained, 
we consider some tests we made to satisfy our- 
selves that the main features of the spectra are not 
introduced by data processing. 

Since the data were obtained by sampling the 
records at 1-min intervals, the shortest period that 
can be included in the spectra is 2 min. Before samp- 
ling the traces, we smoothed them to eliminate any 
periods shorter than 2 min to avoid aliasing. No 
appreciable alteration was involved, however, 
since the amplitudes of components in the 2 min and 
shorter range are small compared with those with 
longer periods. 
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FIGURE 39.—Spectra of the ionospheric disturbance of June 12, 1969, computed with 30-s scalings and 1.0-min scalings, showing that 
the spectra are not appreciably affected by scaling every 1.0 min 


To see whether the spectra were appreciably 
affected by choosing a sampling interval of 1 min, 
we scaled at 30-s intervals the disturbance of June 
12, 1969, recorded at Boulder. The spectrum analy- 
sis was carried out using (1) all scaled points and 
(2) alternate values (1-min intervals). The resulting 
spectra are plotted in figure 39, showing that no 
advantage is gained by making the sampling interval 
less than 1 min. 

Gravity wave components with periods longer 
than about 30 min were also eliminated from the 
data. These are not of interest here and, being of 
relatively large amplitude, tend to obscure the 
acoustic wave components by compressing them 


into the bottoms of the graphs. Hence, the spectra 
are cut off arbitrarily at a period of about 15 min. 

A further test was carried out to see whether 
spurious components are generated by the com- 
puter program. The analytical function 


4 


was used to generate data for a composite signal 
containing 3- and 4-min waves. The resulting spec- 
trum (fig. 40) indicates no spurious spectral com- 
ponents. The line broadening in figure 40 results 
from the finite duration of the function being 
analyzed. 


y=20+ 10 sin (= d +sin (Z ] (35) 
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FIGURE 40.— Spectrum of the function y — 20 + 10{ sin(27/3)t+ 
sin(27/4)t} to test the computer program for spurious 
responses 


Having established that no serious errors are 
introduced in the process of spectrum analysis, 
we now examine the extent that the spectrum 
changes during the disturbance. This is done for 
the disturbance of April 18 in figure 41, which 
shows the spectra for each of the 4 hr from 0500 to 
0900 UT. The evidence suggests that the spectral 
content does change during the disturbance. 

Spectra for the four events described in section 
5 are shown in figures 42, 43, 44, and 45. The 
following features are common to all four: 

1. Most of the acoustic wave energy lies between 
periods of about 2.5 min and about 5 min. This is 
what would be expected for a source in the lower 
atmosphere because the longest period is deter- 
mined by the acoustic cutoff period in the lower 
parts of the atmosphere while the shorter periods 
are suppressed by viscosity (Georges 1969). 

2. There is a range of periods between about 
5 min and about 10 min with relatively little power. 
This range corresponds to the forbidden region of 
figure 7. The shorter period limit is the minimum 
acoustic cutoff period in the atmosphere; and the 
longer limit is the maximum buoyancy period 
that, in this case, is at the level of radio reflection. 
The increase of power for periods greater than 
about 10 min is due to background gravity waves 
nearly always present. 
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FIGURE 41. — Spectral analyses of subintervals for the ionospheric 
disturbance of Apr. 18, 1970, indicating that the spectrum 
changes during the life of the disturbance 


3. The acoustic range is dominated by two power 
peaks, one near 3.5 min and the other near 4.5 
min. The power scale is normalized so that the 
higher of these two peaks is 10. 

The spectral peaks at 3.5 and 4.5 min never 
appear on the Doppler records with sufficient 
strength to be detected by visual inspection when 
there are no thunderstorms within several hundred 
kilometers. To see if these peaks exist but are too 
small to be seen directly in the normal fluctuations 
that occur every day, we analyzed 3 days that were 
free of severe weather (viz, April 11, May 7, and 
June 19, 1970). The peaks were not present on any 
of these days; therefore, there is no evidence that 
their source is anything other than thunderstorms. 
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FIGURE 42. — Spectrum of the ionospheric disturbance of Apr. 18, 1970, at 0500—0900 UT 
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FIGURE 43. — Spectrum of the ionospheric disturbance of Apr. 30, 1970, at 1000-1200 UT 
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FIGURE 44. — Spectrum of the ionospheric disturbance of May 13, 1970, at 0300—0500 UT 
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FIGURE 45. — Spectrum of the ionospheric disturbance of June 11, 1970, at 0300—0700 UT 
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FIGURE 46. —S-shaped Doppler record of the ionospheric disturbance observed on Apr. 25, 1970, at 1600-2030 UT 


Ionosphere 


(0,0) m E | —— (X, O) 


FIGURE 47. —Reflections from a sinusoidal surface 


G. Disturbances Producing S-shaped Records 


The ionospheric disturbance observed on Apr. 
25, 1970, is shown in figure 46; the disturbance is 
quite different from those described in subsections 
5B through 5E. The periodicity of this disturbance 
lies in the range 7 to 10 min, and the waveform 
of the Doppler record is highly distorted in an 
S-shaped fashion. As we have seen in figures 30 and 
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FIGURE 48. — Variation of path length with time 
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FIGURE 49. — Variation of Doppler shift with time 


33, this type of disturbance precedes some of the 
3-min thunderstorm events. We feel that some de- 
scription of this type of disturbance is appropriate. 

The S-shaped records can be explained in terms 
of a corrugated reflecting surface that moves 
horizontally with constant speed. We can make a 
simplified analysis of the problem by postulating 
a sine wave shape to the reflecting surface at a 
mean height H with wavelength L and amplitude 
D (fig. 47). For our analysis, we assume that the 
transmitter and receiver (T and R) coincide and 
move with a speed V while the ionosphere is sta- 
tionary. The former is not quite true in this case, 
but the difference is immaterial. Let the reflecting 
surface be 


z— H +D cos kx (36) 


where k—27/L. When the transmitter is at (X, O) 
(fig. 47), the path (s) up to the point of perpen- 
dicular reflection is given by 


s?= (x-X)*+2 (37) 


where 


0 
X=Vt=7 DUR Dees äi sind ` (38) 


and 0 — kx. Thus, 


s= (H+D cos 0)11 + kD? sin 0}. (39) 


Substituting eq (39) into eq (6) (with P=2s) gives 


di 2fDkV sin 0 


Af WE 40 
C V1+k?D? sin? 0 Si 


for the instantaneous Doppler shift. 

Equations (39) and (40) have been solved nu- 
merically, and typical results for the time variations 
of path length (s) and Doppler shift (Af) are shown 
in figures 48 and 49. We see that, for a certain time 
range and for D > 1 (approximately), there are three 
paths with different lengths and three echoes with 
different frequency shifts. From the S-shaped 
curves, it is possible to get some information about 
the dimensions of the ionospheric disturbance, 
provided that we know the reflection height H 
(e.g., from ionogram analysis) and make several 
approximations. First, let us consider what features 
of the S curves can be measured to a reasonable 
degree of accuracy. These are (1) average period 
T, (2) maximum frequency deviation fmar, and (3) 
the time To from the zero crossing to the turning 
point of the Af(t) curve—see figure 49 (left). 
Since sin 0 ~ 1, we have 
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2fDkV 
A MAX gn 
dE 41) 


from eq (40); and since k?D? <1, 


CT 


di 


Afmax (42) 


where 7(=277/kV ) is the period of the disturbance. 
Equation (42) tells us that, to a first approximation, 
the peak frequency shift is independent of the side- 
way motion of the ray path; therefore, the amplitude 
of the disturbance D can be estimated. A little 
more information can be gleaned from this analysis 
by noting that the approximate condition for a 
switchback to occur is that 


HDE - 1. (43) 


The critical condition exists around the period 
marked A (at 1905 UT) in figure 46. At this time, 
Afmar ~ 0.34 Hz; hence, D=0.69 km. From table 
5, the reflection height H is 179 km that, on substi- 
tution into eq (46), gives k=0.09 km-t, L=70 km, 
and 

Vy = L[r — 166 m/s. (44) 


We can also find Vj from the time displacements 
between the traces in figure 46. In practice, both 
the period and shape of the traces vary from one 
cycle to another; this introduces difficulties in the 
measurement of period and time displacement. The 
zero Doppler frequency crossings, however, are 
essentially independent of the amplitude of the 
reflecting surface. Hence, we use these zero cross- 
ings to determine the average period and average 
time displacements during the approximate time 
interval 1800 to 1845 UT when the S-curves are 
most clearly defined. Cross correlograms of the 
zero crossings from El Reno (E), Shawnee (S), and 
Stillwater (U) are shown in figure 50. Because 
ambiguities occur in matching cycles, we consider 
the following possible time displacements (in s): 


Tes =+ 120 or — 250, 
Tgv =+ 300 or — 100, and 
Tsy =+ 180 or — 240. 


Acceptable combinations of time displacements 


TABLE 5.—Approximate periods and heights of 
disturbances that produced S-shaped records dur- 
ing 1970 


Date Height Period 
(km) (min) 
FT E 179 6.8 
An NE E o oom co 207 9.0 
Té EEN Ne 230 10.5 
must satisfy the relationship 
Tes - Tsv — Ter, (45) 


Applying this criterion to the above data leads to 
the selection of a unique set (viz, Tgs= 120, 
Tsv = 180, and Ter = 300) that also corresponds to 
the highest peaks of the correlograms. Substituting 
these times into eq (12) gives /,-—180 m/s that 
agrees relatively well with eq (44). This selection 
of time displacements leads to an azimuth of 35? 
east of north. As shown in figure 51, this agrees 
excellently with the locations of severe thunder- 
storms. 

Our next step is to estimate the vertical speed 
Vy. From these particular Doppler traces, it was 
impossible to determine Vy for the same time 
interval as for Vy. During the approximate time 
interval 1700 to 1800 ur, the time displacements 
were determined for the different sounding fre- 
quencies. The true reflection heights were obtained 
from the White Sands ionograms. The combined 
data show two possibilities for Vy (viz, 4- 158 m/s 
and — 247 m/s). If, as implied above, the disturbance 
originated in the thunderstorms shown in figure 51, 
phase propagation should be downward (see subsec. 
4B); this leads to the choice of —247 m/s. From 
figure 15, we see that the direction of phase 


propagation is close to the critical value 
$c = {sin (r5/7i) ] where 
dh, —tan^! (Vu/ Vr) = 126°. (46) 


This gives a value of 11.8 min for the intrinsic 
period 7;, corresponding to a buoyancy period of 
9.5 min for the 1962 U.S. Standard Atmosphere. 
Thus, the intrinsic period is considerably greater 
than the observed period, which indicates the 


Kenneth Davies and John E. Jones 35 


APR, 25,1970 1755-1835 UT 
5.1 MHz 


c 


Cross- Correlation Coefficient 


-0.5 


Delay Time, min 


FIGURE 50. —Correlograms of zero crossings of S-shaped records 


presence of a wind component U given by 


dër 
U= e V, — 76 m/s (47) 
i 
(see eq 27). 
An alternative explanation is that there is Vy = 180 m/s 
: 1800 - 1845 UT 
no wind present but that the buoyancy period e = 35° } 


TB =T Sin $.— 5.5 min. Substituting this value into 
eq (14) gives 
H=7.22+25 H' km. (48) 


Near 180 km, we should expect H’ to be in the range 
0.1 to 0.3, so that the scale height should be in the 
range 9.7 km to 15.7 km, corresponding to tempera- 
tures of 310°K to 500°K. Even the larger of these 
temperatures is very much less than the 1280°K 
given by the 1962 U.S. Standard Atmosphere. The 
Standard Atmosphere gives H —36.64 at 180 km, 
which would require that H' — 1.3, an unrealistically 
large gradient. 

Now, the measured speeds of the acoustic dis- 
turbances of section 5 are about 1 km/s and greater, 
suggesting a higher sound speed consistent with 
the U.S. Standard Atmosphere. Hence, we conclude 
that the S-shaped disturbances are evidence that 
the F region has neutral winds with speeds greater FIGURE 51.— Radar summary and horizontal trace velocity for 
than about 70 m/s. Apr. 25, 1970 
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The S-shaped records were observed also on 
two other occasions (table 5), both of which pre- 
ceded 3- to 5-min disturbances. From the three 
events, the indication is that the observed period 
increases with altitude. 


6. POWER IN IONOSPHERIG (NEUTRAL) WAVE 


The power output in the infrasonic range is im- 
portant in identifying the nature of the source. 
With the information available, it is impossible to 
determine the power accurately, but we can get 
an order of magnitude estimate from the arguments 
presented below. 

The energy density in a plane sound wave with 
peak particle velocity v is !/» pv? where p is the 
mean density. Hence, the power flux density is 
1/) pv?u where u is the group velocity. Now we have 
seen that, in the F region, the acoustic cutoff 
period and the buoyancy period are much larger 
than the wave period (x 4.5 min). Therefore, the 
waves are essentially longitudinal, and the group 
speed is very nearly the speed of sound. 

The Doppler technique measures the vertical 
motion of the electrons. This vertical motion is 
related to the neutral particle motion via the geo- 
magnetic field as illustrated in figure 52. The motion 
V of the neutrals is essentially along the k vector, 
but only the component of this motion that is paral- 
lel to the magnetic field is imparted to the elec- 
trons. We observe Vi, the vertical component of 
this electron motion. The angle B between the propa- 
gation vector and the magnetic field is given by 


cos B — sin $ cos] cos Y+ cos sin! (49) 


where J is the dip and w is the azimuth of wave 
propagation with respect to the magnetic meridian. 
The vertical component of electron velocity is 


Vy — V sinl cos 8 (50) 
that causes a Doppler shift: 


aj 2E Ey 2f V sin! cos p. (51) 
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X (Horizontal) 


Magnetic Meridian 


Be Magnetic Field 


FIGURE 52. — Geometry relating motion of the neutral atmosphere 
to the vertical motion of the electrons 


The power density flux is given by 


Ec PIENE SU 
S gee T f sin Í cos 3) (92) 


Let us further assume that we have a point source 
near the ground which radiates isotropically over a 
hemisphere. If our ionospheric observations are 
made at an effective distance r from the source, 
the total radiated power P is then given by 


mr? pCc? Af? 
4f? sin? I (sin $ cos I cos Y+ cos $ sin 1)? 
Po 


- sin? I (sin $ cos I cos Y+ cos $ sin I)? 


(53) 


P=27rrS= 


where Po is the source power that a vertically 
traveling wave must have to produce the same 
Af at the same height where the magnetic field is 
vertical. 

To obtain an idea of the magnitude of Po, we use 
the following atmospheric model: 


p 1071? kg - m^?, 
C — 800 m/s, 

Af — 1 Hz, 
f —5-:105Hz, and 


c—3-10950m«:s^!. 
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FIGURE 53. — Variation of ionospheric response with the direction 
of the ray path for I = 65? 


Since the average height is about 250 km and the 
area of observation lies between horizontal ranges 
of 0 and 200 km from the source, a reasonable 
value of ris 300 km. These values give 


$S9—35-107$W - m? 
and 


The ratio P/Po calculated from eq (53) is near 1.5 
for a vertically traveling acoustic wave (fig. 53) 
and reaches infinity for $—90? and y= 90°. The 
latter condition means that the wave propagates 
horizontally and perpendicularly to the magnetic 
field (magnetic east-west propagation); and under 
these conditions, there will be no electron motion. 
A typical value of P/Po is 2.5, so that the power in 
an acoustic wave source giving rise to an iono- 
spheric disturbance is about 5X106 W. The 
smaller Doppler shifts detected by our Doppler 
system are about 0.1 Hz, which requires a source 
power of about 0.5 x 106 W. 

Figure 53 shows that the ionospheric response is 
sensitive to the position of the source relative to 
the point of observation. The importance of this is 
illustrated in figure 54 for sources in the tropo- 
sphere when the source and observing point lie 
in the magnetic meridian. With the source to the 
south of the ionospheric sounding point, as in 
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FIGURE 54. — Effect of the geomagnetic field on the sensitivity of 
the ionosphere to the location of the source 


figure 54 (left), there will be only a small ionospheric 
disturbance because the acoustic wave tries to move 
the electrons across the magnetic field. With the 
source to the north, as in figure 54 (right), the elec- 
trons will move more easily along the magnetic field. 
Hence, the ionospheric response to a given source 
will be greater when it lies to the north of the sound- 
ing point than when it is to the south. Chang (1969) 
tested this ionospheric sensitivity to source posi- 
tion by using radio data for the path from Long 
Branch, Ill., to Boulder. With the path midpoint as 
the center, Chang divided the circle with 250-km 
radius into a northern semicircle and a southern 
semicircle. For the 14 disturbances he studied, all 
showed thunderstorm activity in the northern 
semicircle; and in 10 of these, the predominant 
activity was in the northern semicircle. In general, 
very few thunderstorms will lie in the magnetic 
meridian and at such locations that the acoustic 
ray path at the radio reflection point is perpendicular 
to the geomagnetic field. Also, the motion of the 
neutral air parcels is slightly elliptical (fig. 8) 
so that there is always some component of motion 
parallel to the geomagnetic field. Hence, very large 
thunderstorms at any location within a radius 
of about 200 km may produce ionospheric 
disturbances. 


7. MICROBAROGRAPH OBSERVATIONS 


If the infrasonic disturbances are indeed pro- 
duced by thunderstorms, we might expect evidence 
of these waves near the ground. To investigate 
this, we installed a microbarograph at the National 
Severe Storms Laboratory in Norman for the period 
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FIGURE 55.—Sample microbarograph records (Norman, Okla., on June 11, 1970, at 0300-0700 UT) showing scaling procedure 
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FIGURE 56. —Microbarograph spectra 


April-July 1970. The microbarograph was loaned 
by Dr. H. Montes of Teledyne Isotopes, Westwood, 
NJ. 

Sample microbarograph records are shown in 
figure 55. Several records have been spectrum 
analyzed by the technique described in section 6, 
and the spectra are presented in figure 56. In gen- 
eral, we do not find convincing evidence for the 
presence of strong signals with periods near 3.5 
and 4.5 min. This does not mean that signals with 
these periods were not present. It is quite possible 
that they were swamped by the background. Much 
of this background is filtered out of the ionospheric 
disturbance by the propagation characteristics of 
the atmosphere. 

Evidence that infrasonic waves originate in some 
of the thunderstorms associated with ionospheric 
disturbances has been provided by Vernon Goerke 
of NOAA. Using spaced microphones at Boulder, 
Goerke and Woodward (1966) monitored the direc- 
tion of arrival of infrasound with periods between 
20 and 50 s. Figure 57 shows radar summaries coin- 
ciding with two ionospheric events in 1969 and the 
azimuths of arrival of infrasound at Boulder. Hence, 
despite the negative results in the microbarograph 
recordings, there is reason to believe that infra- 
sound is generated by thunderstorms that produce 
ionospheric disturbances. 
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FIGURE 57.— Radar summaries and azimuths of infrasound received at Boulder (con. on p. 40). 
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8. SOME STATISTICAL ASPECTS 


Tables 6 and 7 list all the acoustic wave events 
detected at Norman, Okla., during 1969 and 1970. 
The disturbances have been designated as minor 
and major, depending on the magnitude of the rms 
Doppler shift— minor < 0.1 Hz, major > 0.1 Hz. 
The tables show that most disturbances are near 
the limit of detectability of our Doppler system. 
Some of the major events may start and end as 
minor events. Other major events appear to start 
rather abruptly. The sferics were recorded at 
Norman through April and May of 1970. From the 
daily geomagnetic and solar data listed in tables 6 


and 7, we see little or no evidence of correlation with 
the occurrence of the acoustic disturbances. An ex- 
amination of the Kp values for the hours covering 
the observed disturbances and 5 quiet days and 5 
disturbed days of the months have likewise failed to 
yield any sign of geomagnetic influence. 

We have also sought to find temporal variations in 
the characteristics of the spectra. For example in 
figure 58, we have plotted the frequencies with 
maximum powers as functions of season and uni- 
versal time. There is no convincing evidence of 
any systematic change in these frequencies. In 
figure 58 (bottom), notice that the lower frequency 
points show little dispersion. This is confined by 
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TABLE 6.— Acoustic wave events during May 9— TABLE 7.— Acoustic wave events during Mar. 4— 
Sept. 25, 1969 (f — 3.3 MHz) July 30, 1970 (f— 5.1 MHz) 
Minor Major Minor Major 
D Mag- 
Date Begin End Begin End Magnetic — 2800-MHz Date Begin End Begin End Sferics netic 2800-MHz 
(UT) (ur) (UT) (UT) (Ap) flux (UT) (UT) (UT) (UT) (Ap) flux 
May 15. 02 10 131 159 Apr. 
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June 1 07 09 5 113 18 0 16 05 0 * ` I8 142 
June 12 00 08 03 06 18 237 19 ll 16 * 2] 136 
Junel3 04 07 15 229 ot. 30 28 90 126 
June 18 02 08 4 154 23 15 20 * 14 128 
June 20 10 13 10 145 24 fa 12 16 130 
June 22 o 2 133 16 18 
June 23 00 03 5 132 91 La 12 * 17 135 
June 24 00 12 0 06 12 129 03 05 10 138 
June 25 ll 14 10 120 eg 02 ji i 3s : 18 ES 
lue 97? "y 08 5 114 May 
June 29 00 04 00 0 4 118 50 MER A E ee 
No records July 18- August 15 1] Kë 03 3 176 
Aug. 27 12 14 21 168 SP 
12 09 12 15 180 
03 12 03 05 6 193 
13 03 12 09 11 
22 23 
15 02 08 10 206 
05 09 45 154 
- | 14 16 
55 : 29 ll 14 * 13 159 
30 00 04 à H 164 
9.0 June 
45 10 [po 12 6 146 
22 07 03 07 
40 H 2. E UD . 5 172 
35 13 04 06 04 06 10 200 
a i 14 06 12 8 207 
a OL ir 2 a 0 X 3 10 20 3 9 19 29 ee NER ^ 
A-| March F wi "Was "| "June | at M ae ls x 
2 eg (9 04 00 04 6 149 
c 08 12 
er 91 þe 08 35 160 
k jà- ud 
July 
2 14 16 11 189 
11 14 16 10 138 
12 01 05 14 135 
15 d 07 01 05 9 122 
05 07 05 07 


Ww 02 mw —o7.0 IFR 
UT 


* Events that coincided with high rates of electromagnetic 
pulsations caused by lightning discharges from thunderstorms 
FIGURE 58. — Seasonal and diurnal variations of the power spectra within a radius of several hundred kilometers of Norman, Okla. 

peaks (Taylor 1971) 
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FIGURE 59.—Histogram showing the distributions of the peaks 
in the spectra of ionospheric disturbances 


the histogram in figure 59 showing that most of the 
observed values are at 3.65 mHz (r ~ 4.6 min); 
there are no events with smaller values. On the 
other hand, the high-frequency peak (~ 4.7 mHz) 
is more variable. These facts reflect the general 
properties of the spectra (viz, that the peak near 
4.6 min is narrow, whereas the peak near 3.5 min 
is broad). There seems to be no fixed relationship 
between the two frequencies. 


9. DISCUSSION 


The experiments carried out in Oklahoma prove 
that ionospheric disturbances with periods in the 
range of 2 to 5 min are indeed generated by thunder- 
storms. In fact, during the thunderstorm season, 
the ionospheric radio data and the radar summaries 
were often collected independently; and by inspect- 
ing one set of data, we could predict the properties 
of the other. 

The interesting questions to ask are (1) “How 
is the acoustic wave produced?" and (2) “How does 
it propagate from the source to the ionosphere?” 
In the treatment here, we have assumed that the 
disturbance is generated by the motion of the 
thunderstorm and propagates to the ionosphere 
as an acoustic wave. However, other possibilities 
should be considered, for example, (1) the dis- 
turbance could arise from a shock wave caused 
by the heating effect of lightning and filtered by 
dispersive propagation or (2) Herman (1970) has 
suggested that high frequency radio noise pulses 


(sferics) emitted by lightning discharges propagate 
to the F region where they are absorbed. Utlaut 
(1970) and his collaborators have shown that a 
certain type of ionospheric disturbance (spread F) 
can be produced by heating the F layer with a 
powerful radio transmitter. 

Before discussing the various possibilities, let 
us summarize the facts that any theory must 
explain. 

1. Disturbances occur only when severe thunder- 
storms are in the vicinity (~ 200 km). 

2. Not all thunderstorms with tops higher than 
40,000 ft and within a radius of 200 km produce 
detectable disturbances. 

3. The spectra of the disturbances have two peaks, 
one (near 4.5 min) that is narrow and frequency 
stable and the other (near 3.5 min) that is broader 
than the first and tends to vary somewhat in fre- 
quency from one event to another. 

4. The direction of propagation of the ionospheric 
disturbance may change abruptly during a long- 
lived event. 

A likely source of these waves lies in the me- 
chanical motions of the air mass that constitutes a 
thunderstorm cell. Based on photographs of clouds, 
evidence for oscillatory motions in thunderclouds 
over Arizona has been advanced by Anderson 
(1960). His data show a dominant period of nearly 
10 min. Superimposed on the low frequency are 
numerous higher frequency oscillations with periods 
near 2 min. If anything, Anderson’s spectra indicate 
a minimum of energy in the 3- to 4-min range, 
although most of the cloud tops studied by Ander- 
son are considerably lower than the 40,000-ft 
requirement established by Baker and Davies (1969). 
Furthermore, in contrast to Anderson’s results, 
our power spectra have no peaks near 10 min. A 
period near 10 min is consistent with buoyancy 
oscillations in a moist atmosphere. We are, there- 
fore, led to believe that the disturbances with 
periods near 3.5 and 4.5 min are not produced by 
simple buoyancy oscillations. 

Our data suggest that at least part of the wave 
spectrum is produced by elastic oscillations at 
the tropopause level. These may be produced in 
the storm cells themselves, or the cells may excite 
natural oscillations in the surrounding gas at the 
tropopause. 

The temperature gradient at the tropopause is 
zero; hence, the acoustic cutoff period is related to 


Kenneth Davies and John E. Jones 43 


the temperature by 


Aart. 


Yg 


=18.39 VT (54) 


Ta = 


where Ta is in seconds and T is in degrees Kelvin. 
The mean temperature over Oklahoma in 1969 at 
the 200-mb level, which is approximately the tropo- 
pause height, was — 55°C or 218°K. This corresponds 
to an acoustic cutoff period of 4.52 min, which is 
very close to one of the peaks in the observed 
spectra. The sensitivity of Ta to the temperature 
fluctuation AT is given by 


ATa=— 9.2 T1? AT. (55) 


The 200-mb temperature extremes over Oklahoma 
are — 55° 10°C so that Ar=+0.1 min. This would 
explain the observed stability in the 4.5-min peak. 
At this point, it is of interest to calculate the 
pressure perturbation (p'/po) in the ionosphere: 
D yV y 


po cr” 


(56) 


where u is the horizontal velocity of the air parcels. 
Taking Vy=1 km/s, C— 0.8 km/s, and u— 10 m/s, 
we find that 


-= 0.022 = 2.2%. (57) 


0 


Remember that the power in the acoustic wave is 
only one-millionth of the power generated by a big 


storm. 
Another possible cause of these disturbances lies 


in the shock waves generated by lightning dis- 
charges (Few 1969). Estimates of the power gen- 
erated by lightning in a severe thunderstorm vary 
over a wide range from about 109 to about 10!'? kW 
(Loeb 1954 and Vonnegut 1960). Taking a value of 
10? kW, we see that there is adequate power avail- 
able even when allowing for a highly inefficient 
conversion process. Nevertheless, it is difficult to 
see how such a process can produce the spectrum 
observed in the ionosphere. 

Our data suggest that the sources of the iono- 
spheric disturbances are the individual thunder- 
storm cells rather than the storm system as a 
whole. Evidence for this is provided by the sudden 
change in direction of the horizontal trace velocity 


on April 30. It appears that the waves from one 
source die away while those from another source 
build up and eventually dominate. Further support 
for this hypothesis comes from the ground-based 
infrasonic recordings. These observations of distant 
thunderstorms show rather abrupt changes of azi- 
muth indicating changes of sources. This generally 
agrees with the meteorological classification that 
divides the life of a cell into a buildup stage, a 
mature stage, and a decay stage (Palmén and 
Newton 1969, ch. 13). The average life of such cells 
is about 0.5 hr with large “‘super cells" persisting 
for several hours. These lifetimes are consistent 
with the durations of ionospheric disturbances. 


10. CONTINUING OBSERVATIONS 


During 1971, the recording of ionospheric dis- 
turbances has continued with observations made in 
Florida. A single 5.2-MHz transmitter was set up 
at a National Weather Service field site near Lake- 
land, with the receiver located at the University of 
Florida campus in Gainesville. We intended to 
look for ionospheric disturbances associated with 
hurricanes and to see if coastal thunderstorms 
produce effects that are different from those pro- 
duced inland, such as in Oklahoma. 

Records obtained from April through August 1971 
do not show nearly as much acoustic wave activity 
as records of 1970 in Oklahoma. In fact, there was 
not even one ionospheric disturbance that would 
have been classified as a major event, or one that 
had a long-enough duration for spectrum analysis. 
Although the data have not yet been thoroughly 
analyzed, the difference may be because the 
tropopause over Florida is higher (about 16 km) 
than that over Oklahoma (about 12 km). Therefore, 
if penetration of the tropopause is necessary to 
produce disturbances, higher cloud tops would be 
required in Florida; and the 40,000-ft criterion 
would not apply. 
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